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The probabilitiy distribution density of random values of squared 
functional on Wiener process trajectories 

Virchenko Yu.P., Vitokhina N.N. 

1. We consider the problem about the distribution density calculation of 
the random variable Jt[w] in this work. Here 

T 

Jt[x] = I x\t)dt (1) 



o 

O; is the functional in the L2([0,r]) space and {w{t);t G [0,r]}, T > are 

; trajectories of the standard Wiener process on [0,r], i.e. Ew'^{t) = t. This 

problem is classical as well as analogous problems concerns the probability 
distribution density calculation of random variables performed by additive 
^ ; squared functionals on trajectories of gaussian random processes. For such 

i processes, the calculation problem connected with characteristic functions of 

lO ! random values under consideration is solved principally. Up to present time, 

■ a great number specific problems are solved on the basis of this method which 

Ph; have various applications (see [1], [3]). For example, such a result concerns 

^ ; the normal markovian process (the Ornstein-Uhlenbeck process) has been 



still obtained in the work [2]. However, the problem of restoration of the 
distribution density on the basis of obtained characteristic functions remains 
weakly investigated in the sense of the constructing of approximations with 
^ ; the guaranteed accuracy being suitable for use whereas in the area of the 

random variable changing. The usual approach to the solution of this problem 
(see, for example, [4]) results in some approximated formulas for distribution 
densities suitable for the estimation of great fluctuation probabilities, i.e. 
in the asymptotic area x ^ oo of the random value changing. Here, we 
study the calculation problem of successive approximations of the probability 
distribution density f{x) connected with random values of the functional (1) 
defined in any compact interval [0, M], M > 0. 

2. Since trajectories {w{t);t > 0} of the standard Wiener process Ew^{t) = 
t are continuous with the probability one, then the random value Jt[w] is de- 
termined almost sure for each of them. 

The generating function of the random value (1) is evaluated by the for- 
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mula (see, for example, [1]) 
QriX) = Eexp(-AJTH) = 



ch 



-1/2 



Qt(A) = QiiXT') 



The distribution density f{x) of the random value Jt[w] is defined by the 
inverse Laplace transformation 
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ch f X'/'T 



-1/2 



dX 



(2) 



where c >0 and, for the integrated function, the cross-cut in the complex 
plane A is done along the negative part of the real axe. 

It is convenient to introduce the density g{x) = T'^ f\T'^x). The replace- 
ment of the integration variable in Eq. (2) gives 



9{x) = 
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ioo+c 



27ri 



ch ( Ai/2 



-1/2 



dX = 
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ioo+c 
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exp(2Aa: - AV2) 
7!^ J Vl + exp(-2A-i/2) 

—ioo+c 



1/2 



dX 



(3) 



in this case. Let us prove the following theorem. 

Theorem. The density g{x) is represented by the following absolutely 
converging series 
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(4) 



the Nth remainder of this series is estimated by the value 



\9ix) - g„.,{x)\ <J^ (iS^) + 1/4) exp (- + 



TTX^ 



X 



(5) 



□ We shall put c = in Eq.(3) since singularities of integrated function 
are on the negative part of the real axe. We deform the integration contour 
to the contour C consisting of the consecutive transitions of following ways 
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{s -i£;s e (-00; 0]}, {se"-'; s E [-7r/2; 7r/2]}, {-s + is; s G [0; +00)}. Such 
a deformation is permissible since 

ch (aV2) I' = ch (a^/^) ch ((A*)^/2) = i [ch (2Re(A^/2)) + ch (2ilmiX'/') 

> i (ch (2R^/'^cos{(p/2)^ ~ ^) = (i?^/^cos((^/2) 

where A = i^e*"^ and, on the arch {A; ip G ['7r/2; tt)} of the circle, the following 
estimation of the integrated expression module in Eq. (3) is valid 



exp(Ax) 



(ch(Ai/2)) 



1/2 



< 



exp(a:i? cos ip) 

[sh(i?i/2cos((^/2))] 



1/2 



It guarantees the fulfillment of the Jordan condition at a: > on the arch 
i?^/^ cos((/?/2) < £ at any small e: > since cos (p < 0. The same takes place 
for the arch {A; (p G (— 7r,7r/2]}. 

In Eq.(3) we shall realize the replacement of the integration variable A^/^ = 
q, then A = q^, dX = 2qdq. Thus, the contour C in the complex plane A will 
be transformed to the line and, after the transition to the limit £ ^ 0, it will 
be transformed to the line {q = is; s G M} in the complex plane q. After 
these transformations, we have 



> 



g{x) = 



V2 



too 



q 



exp(2g ,7' — q) 
, 1 + exp(-2g) 



1/2 



dq 



—too 



Now, we pass to the integration on the variable s,q = c-\-is,dq = ids. Then 
we obtain 



V2 f^ eM-2o^s^-^s) Y\^_ 
^ 7T J \ l + exp(-2zs) J ^ ^ 

> s of the integra- 



In the last integral, we shall produce the shift s + z(4x) ^ 
tion variable. Therefore, we obtain 



g{x) = i 



.V2 



+00 



TT 



exp(-2x(s - z(4x)-i)2 - (8x)-i) 
1 + exp(-2z(s + z(4x)-i) - (2x)-i; 



1/2 



ds = 



—00 
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+00 



exp(— 2xs^) 



1/2 



+ e-V(2a;) exp(-2is) 



ds. (7) 
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We decompose the denominator of the integrated expression in Eq.(7) into 
the series converged at any x > and at any s G M 
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(1 + e-^/^2^^exp(-2is))-i/2 = ^(-1)^^?^— exp(-//(2x)) exp(-2i/s). 



l=Q 



The convergence is uniform in any area [0, M] x M in the plane (x, s), M > 0. 
Substituting the last expression in Eq.(7), we obtain 



-\/2 f 
g{x) = — exp(— (16a;)~^) / (is + (4a;)~"'') exp(— xs^) x 
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exp(— 2z/s) exp(— //(2x)) 



ds = 
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exp(— (16x) ^) 



(20! 



+00 



X exp {-x{s + iZ/x)^ - + 1/2))] ds . 

The transposition of summation and integration is based on the uniform 
convergence of the series on s at any fixed x. 

In each summand of the sum, we shall produce the shift s + il/x ^ s of 
the integration variable, 

g{x) = — exp(— (16x) ^) x 

TT 
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{is + x-'^ [I + 1/4]) exp (-xs^ - x'Hil + 1/2)) 



ds 
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Let us represent the integral as the sum of two integrals in accordance with 
the expression in the bracket before exponent. The integral corresponding 
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the summand is converts to the zero 

+00 

2 



i J |;(-l)'^exp(-xs^)exp(-x-iJ(J + l/2)) 

-OO 



ds = , 



due to the oddness of the integrand function. The integral corresponding the 
summand [/ + 1/4], transforms as fohows 
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+00 r- ^ 

J exipi-xs')ds J2^-iy^2 + V4] exp + 1/2)) 

-oo _ '- 



Substitution of this expression in Eq. (8) taking into account the value of the 
Poisson integral results in the formula (4). 

Since the series (4) is alternating in sign then the remainder of the se- 
ries does not exceed the first summand among rejected ones. Hence, the 
estimation (5) is valid. ■ 

Corollary. The following estimation takes place 



(9) 



□ Let us estimate the remainder of the series (4). For this purpose, on the 
basis Eq.(4), we write down the density g{x) in the form 

[2 °° 



1=0 



where 



ai 



(20! 
4^(/!)2 



(/ + 1/4), hi{x) = x~^/^ exp 



(/ + 1/4)^ 
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Also, we may find maximums on x of the functions h]s[{x), n = 1,2,3, ... . 
Equating to zero the derivative on x of this function 

h'^ix) =x~'^hN{x) [(iV + 1/4)2 -3a:/2] 

we find the solution of this equation. It is the point of the maximum of 
the function h^ix) being unique for each N ^ 



= 3(^+1/4)^ 
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Hence, the estimation of the Nth remainder is 



, / N / M /2 , , ^ 3 / 3 / (2iV)! 

\9{x) -gN-i{x)\ < \ -aNhN{x^) = 7;\ — 

\ IT Z \ ire 



V4^(iV!)^ 



{N + 1/4) 
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Further, we estimate the coefficient having done more transparent the 
obtained estimation. It is made by the following way 



= exp 



N 



(2Ar- 1)!! 



^ln(l-(2/)-i) 



1=1 
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n 

1=1 



< exp 



•2/-1' 
21 



= n(-i) = 



2^ 

-1/2 







< exp 


-i(l + ln TV) 



N 

in view of validity of inequalities ln(l — x) < —x at > and 



^ 1 
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> 1 + 



dx 
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= 1 + In 



Since y^S/vr < 1 then Eq.(9) takes place. ■ 

3. Now, we estimate the approximation accuracy of probabilities 
PrjJyfw] > c} which are obtained on the basis of functions gN{-)i N = 
1, 2, ... . Since f{x) = T-^g{T-^x) then 

c 

Ft{Jt[w] >c} = J g {T-^x) dx = 1 - R{c) , 



where 



c/T2 

i?iv(c) = J gN-i{x)dx 
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Designating the righthand side of the inequahty (5) by Qn{x)i we have 
\g{x) — gN-i{x)\ < Qn{x). Further, we determine the function 

c/r2 

Pn{c) = 1 - Rn{c) , Rn{c) = J gN-i{x)dx. 



Our problem is the reception of the top estimation of the deviation 
|Pr{J7-[w;] > c} — P/v(c)|. From the inequahty (5), it follows 

-Qn[x) < g{x) - gN-i{x) < Qn{x) . 

Integrating between limits and c/T^, we obtain 

c/T2 c/T2 c/T2 

— J QN{x)dx < R{c) — J gN-i{x)dx < J Q]sf{x)dx . 



Hence, 



\Ric)-RNic)\ = 



c/T^ 



{g[x) - gN-i{x))dx 



< j QN{x)dx 



This gives the desired estimation 

\Pt{Jt[w]>c}-Pn{c)\ = \R{c)-Rn{c)\ < j QN[x)dx 



(10) 



At last, we calculate the integral in the righthand side. Due to this, the 
estimation (10) becomes obvious 



j QN{x)dx = aN\j^ j exp ^ 



[N + 1/Af\ dx 



X ) 



Replacement of the integration variable y = dy = —dx/2x'^/'^ results in 

the formula 



QN{x)dx = 



8 CLjSf 

7r7VTT/4 
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r(Jv+i/4) 

cV2 



■r(A^ + i/4) 



From here, using Erfc(a:) < 1, we find the following estimation being uniform 
on parameters c and T 

/ QNix)dx < y^- 



More exact estimation which takes into account the order of parameters c and 
T is obtained by using the standard inequality Erfc(x) < exp(— x^), 

J QN{x)dx < 



< (™V^)"'exp(-(TiV)Vo). 
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